Fluctuation interactions of colloidal particles 



T. Ocampo-Delgado^ and B. Ivlev^'^ 
^Instituto de Fisica, Universidad Autonoma de San Luis Potosi San Luis Potosi, 
San Luis Potest 78000 Mexico 
and 

^Department of Physics and Astronomy and NanoCenter, 
University of South Carolina, 
Columbia, South Carolina 29208, USA 



For like charged colloidal particles two mechanisms of attraction between them survive when the 
interparticle distance is larger than the Debye screening length. One of them is the conventional 
van der Waals attraction and the second one is the attraction mechanism mediated by thermal 
fluctuations of particle position. The latter is related to an effective variable mass (Euler mass) 
of the particles produced by a fluid motion involved. The most stronger attraction potential (up 
to the value of temperature T) corresponds to a situation of uncharged particles and a relatively 
large Debye screening length. In this case the third mechanism of attraction enters the game, ft is 
mediated by thermal fluctuations of fluid density. 

PACS numbers: 05.20.Jj, 82.70.Dd 



I. THE SHORT STORY 

Systems of charged colloidal particles exhibit a vari- 
ety of unusual physical properties [H, 13, Q- Colloidal 
particles can be arranged into crystal 14] and into struc- 
tures with clusters and voids 0, [1, Colloidal sys- 
tems may undergo different types of phases transition 
d, d, [13, [m, [13, Ell- Topological phase transitions in 
two-dimensional system of colloidal particles have been 
discussed in [3, • Unusual ensembles of colloidal par- 
ticles were observed in Refs. [ll, [13, [l3| . In [Hi buck- 
ling instabilities in confined colloidal crystal were ana- 
lyzed. Interesting behaviors of colloidal in external fields 
were reported in [1^. Colloidal particles accept in an 
electrolyte some surface charge, screened by counteri- 
ons at Debye's length Ad, which results in the repulsion 
potential of Derjaguin, Landau, Verwey, and Overbeek 
(DLVO) [H, [2. The DLVO theory, as a result of solu- 
tion of linearized Poisson-Boltzmann equation, has been 
questioned in [2l|, [23] . The generalization of DLVO in- 
teraction via a modification of counterion screening was 
reported in [2^. 

The long range attraction of like charged particles is a 
matter of challenge and controversy in colloidal science. 
Very schematically, the story of the long range attraction 
is the following. 

Initially, the long range attraction of colloidal parti- 
cles on the micron scale has been reported in Ref. [2^ . 
The authors reconstructed a pair potential from mea- 
sured correlation functions of a large ensemble of colloidal 
particles. 

The authors of Refs. [13, [13, [13| published even more 
surprising results when the attraction was extended over 
a distance of almost 4 /im. They used a laser technique 
releasing two particles and observing their behavior. It 
happened that their measurements were misinterpreted. 
It was observed just a macroscopic hydromechanical ef- 
fect associated with specificity of measurements but not 



a microscopic mechanism of attraction ^28*]. 

Bechinger and his group shown for a high concentrated 
particle ensemble that different types of interparticle po- 
tential result in almost identical correlation functions 
29, 30j. The conclusion was that an interaction potential 
should be drawn from experiments with a pair of colloidal 
particles but not from a large ensemble of them. Another 
difhculty is connected to an uncertainty in observed par- 
ticle position caused by diffraction which can result in 
errors in calculated pair potential. See Ref. [Slj and ref- 
erences therein. To reduce the diffraction uncertainty 
one should use ultraviolet observations. So reliability of 
calculation of pair potential on the basis of statistical 
properties of large ensembles of particles is questionable. 

Despite a lack of reliability of various experiments on 
determination of a microscopic mechanism of attraction 
of like charges it is always intriguing whether they attract 
or not in reality. We consider the following microscopic 
mechanism of attraction of two colloidal particles or a 
particle and a wall. 

1. The conventional van der Waals attraction Uydw 
[131 mediated by electromagnetic fluctuations. 

2. The attraction Ucom, mediated by thermal compres- 
sion fluctuations of a fluid. 

3. The attraction / associated with thermal fluctu- 
ations of particle positions in a fluid. It results from 
variable particle masses (Euler masses) depending on a 
distance between them. This mechanism and the term 
"Euler mass" was proposed in Ref. fs^. 

In the first mechanism an energy of fluctuating electro- 
magnetic waves depends on a distance between particles 
and, therefore, it leads to a force. 

The second mechanism is generic with the first one 
since electromagnetic waves are just substituted by hy- 
drodynamic ones. 

In the third mechanism moving particles drag a part 
of a fluid. A mass of the involved fluid depends on a dis- 
tance between the particles. There is a thermal drift of 
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the particles into the region with a larger effective mass 
which is analogous to classical mechanics. This can be in- 
terpreted as an effective interaction mediated by thermal 
fluctuations of particle positions in a fluid. In formation 
of the variable mass high frequency fluctuations of the 
particles are involved when dissipative hydrodynamic ef- 
fects are not important [sj. This corresponds to Euler 
hydrodynamics and gives rise to the term "Euler mass" . 
The formation of Euler mass reminds the equipartition 
law when a mean kinetic energy is T/2 regardless of dis- 
sipation. This is due to that the mean kinetic energy is 
also determined by high frequencies. 

There is a substantial difference among the above at- 
traction mechanisms. The first and the third ones survive 
when an interparticle distance becomes larger than the 
Debye screening length. On this distance the DLVO re- 
pulsion is very small and the above two mechanisms are 
the only interaction. The second mechanism works when 
Coulomb effects in the fluid are not pronounced. Namely, 
the particles are not charged and the Debye screening 
length is larger than the interparticle distance. A rela- 
tive role of attraction I and Uydw was also analyzed in 
Ref. [H. 

A goal of this paper is to study the three above attrac- 
tion potentials for two particles and for one near a wall 
(walls). 



II. VAN DER WAALS INTERACTION 

The energy of fluctuating electromagnetic field around 
two particles depends on a distance R between them and, 
therefore, results in an interaction force. It is called van 
der Waals force [13]. This force is mainly determined 
by a typical wave length A of the fluctuating electromag- 
netic field since it should be A ~ i?. The permittivity 
of the particles material s{uj) depends on the typically 
frequency u ~ c/R. We consider relatively large inter- 
particle distance 



dielectric constant Eq for a surrounded medium. When 
the two dielectric constants are close to each other 
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one can use the approach of pairwise summation to calcu- 
late the energy of electromagnetic fluctuations (the van 
der Waals interaction energy) [32, 36] 
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In Eq. (Ill) the integrations occur inside the volumes of 
two bodies. For two identical spherical particles of the 
radius a and center-to-center distance R the integration 
in Eq. g]) results in [H HI] 
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For a spherical particle near a fiat infinite wall (the 
both with the dielectric constant e) Eq. (|4]) yields 
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In Eq. ^ h is the center-to-wall distance. 

In the region of visible light for water Eq ~ 1.77 and 
for usually used polystyrene colloidal particles e ~ 2.40. 
So the parameter (e — eo)/£o — 0.35 can be considered 
as relatively small and, therefore, the equations ^ and 
are reasonable approximations for van der Waals in- 
teraction. 



III. INTERACTION MEDIATED BY 
COMPRESSION FLUCTUATIONS OF THE 
FLUID 
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where ujq corresponds to an absorption peak of s{uj). For 
example, for water ujq ^ 10^^ sec^^. The estimate ([T]) 
reads 100 A < R. 

At a finite temperature the typical wave length of the 
fiuctuating electromagnetic field is hc/T. Below the in- 
terparticle distance is not too large 
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which is equivalent to i? < 7.4 fim at room temperature. 

We consider a typical interparticle distance R of the 
order of one or two microns which relates to the condi- 
tions ll]) and In the optical interval of uj the per- 
mittivity is determined by a refractive index and can be 
substituted by a dielectric constant e for particles and a 



Suppose two particles to be totally fixed inside a hy- 
drodynamic medium and they serve only as obstacles for 
a fluid motion. There is no macroscopic motion in the 
system and the only motion is caused by thermal fluctu- 
ations of the fluid velocity v{r,t). In this case the free 
energy of thermal fluctuations of the fluid F{R) depends 
on the distance R between the particles. The function 



Uco,niR) - F{R) ~ F(c») 
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is an interaction mediated by compression fiuctuations of 
the fiuid analogously to the conventional van der Waals 
interaction mediated by electromagnetic fiuctuations. To 
find the free energy of thermal fiuctuations of the fiuid 
one can start with the linearized Navier-Stokes equation 

M 
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There are two types of fluid motion, one of them is a 
transverse diffusion and the second one is longitudinal 
sound waves associated with the density variation. The 
equilibrium free energy of transverse motions is deter- 
mined by the Boltzmann distribution of their kinetic en- 
ergies and does not depend on the friction coefficient in 
the thermal limit. The free energy of transverse motions 
in the thermal limit depends on the total volume, but 
not on relative positions of bodies. Therefore, transverse 
fluctuations do not result in an interaction. 

Quite opposite situation occurs for longitudinal mo- 
tions, when the total free energy is a sum of energies 
of different sound modes. The spectrum of sound waves 
depends on the distance between bodies R due to hy- 
drodynamic boundary conditions on body surfaces and 
this results in i?-dependence of the free energy. Hence, 
the fluctuation interaction between bodies is mediated 
by hydrodynamic sound waves like the conventional van 
der Waals interaction is mediated by fluctuations of elec- 
tromagnetic ones. Putting v = \Jd4>/dt, one can obtain 
from Eq. (O 



5^' 
9*2 



(9) 



Through thermodynamic relations and the continuity 
equation one can obtain 5p = —ps^'^'^cf), where s is the 
adiabatic sound velocity [3J]. At the typical frequency 
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la ~ 



1 /xm is a particle radius), in- 



volved into the problem, the dissipative term in Eq. Q 
is small and one can write 
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According to the limit of a small friction, the boundary 
condition for the normal derivative = to Eq. (fTO|) 
corresponds to the Euler fluid [13] ■ From general point of 
view, the free energy of a system of harmonic oscillators 
does not depend on friction in the thermal limit. 

In an electrolyte the dispersion law of sound waves can 
be approximated as 
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where Xjj is the Debye screening length. Let us consider 
flrst the case of two inflnite parallel walls separated by 
the distance R. The free energy per unit area of the 
system is expressed as a sum of energies of independent 
oscillators according to general rules of statistical physics 
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Performing the same steps as in Ref. jSJ] , we obtain the 
interaction mediated by compression fluctuations in the 
form 



1 — exp 



-^z + ARyX 



(13) 



In the limiting cases Eq. (|T3|) reads 
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where C(3) — 1-202 is the Riemann zeta function. 
At a large R the interaction ([Tl)) is screened on the 
length Ad/2. A possibility of interaction, mediated by 
non-electromagnetic fluctuations, has been proposed by 
Dzyaloshinskii, Lifshitz, and Pitaevskii [32]. The flrst 
formula (|14p was obtained in Ref. [s^. It is similar to 
the result of Ref. ^36*] for electromagnetic fluctuations and 
perfectly conducting planes. 

The fluctuation interaction (fT3|) depends on the Debye 
screening length. This happens due to summation over 
all wave vectors in the free energy. This results in its 
dependence on density of states which, in turn, depends 
on a form of the spectrum. In our case the spectrum 
u; = s{k^ + A^^)-^/2 brings the Ajj-dependence into the 
free free energy. 

When two objects are not flat but they are close 
enough interacting by small parts of their surfaces which 
are almost flat, one can derive an interaction potential 
from the flat approximation l|14p integrating over the sur- 
faces. For example, for a particle close to a flat wall when 
the center-to-wall distance h is smaller than the particle 
radius a, the interaction can be calculated as in Ref. [ssj 
and has the form 
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The analogous result for two spheres with the center-to- 
center distance R is 



C(3)y_ a 



16 R-2a 



(i?-2a)<a,AD. (16) 



Equations and (fTB)) hold in a non-electrolytic fluid or 
in an electrolyte with a sufficiently large screening length. 



IV. WHICH MECHANISM SURVIVES IN AN 
ELECTROLYTE 

In an electrolyte the interaction Ucom, mediated by 
plasmons, is strongly screened (fT4|) . The van der Waals 
interaction Uydw, rnediated by photons of a frequency in 
the visible range, is not sensitive to the plasmon effects. 
The mechanism of variable mass is connected solely with 
incompressible fluid fluctuations. Therefore, only the 
conventional van der Waals interaction Uvdw and vari- 
able mass mechanisms / can survive in an electrolyte. 
The interaction due to variable mass is considered in the 
following sections. 



V. VARIABLE MASS MECHANISM 

For the purpose of illustration of the variable mass 
mechanism one can consider a simple mechanical analogy. 
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Suppose a classical non-dissipative particle of the total 
energy E moves in the harmonic potential ax^. When 
the particle mass is a constant the mean displacement 
(x) = in the harmonic potential. But in the case of 
a variable mass m{x) the mean displacement (x) ^ 0. 
This is due to that a particle velocity is smaller at a 
region of a larger mass just to keep the total energy to 
be constant. According to that, the particle spends more 
time at the region of a larger mass. This is equivalent to 
the certain effective attraction I{x) to a region of a larger 
mass. As shown in Ref. [s^, for a slow varying m{x) the 
total effective potential becomes ax^ — {E/2) lnm(a;). 

A real particle with friction participates the Brown- 
ian motion characterized by the certain temperature T. 
We briefly repeat here the main arguments leading to 
an effective potential / [3^ . The Langevin equation, de- 
scribing such processes, has the form 
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■qx = stochastic force. (17) 



Short time fluctuations of the velocity x are well sepa- 
rated from the slow drift in an effective potential. In- 
deed, according to fluctuation dissipation theorem, the 
mean value of the kinetic energy (mi^/2) = T/2 cor- 
responds to the equipartition law and is contributed by 
short time fluctuations related to the infinitely large cir- 
cle in the plane of complex frequency. Substitution of 
that mean value to Eq. (fT7|) leads to the effective poten- 
tial V{x) -I- I{x) [H, where 
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The expression (llSp is an exact result in the thermal limit 
(no quantum fluctuations) as the equipartition law. 

We remind what happens in the conventional case of 
position-independent masses. In that situation in the 
thermal limit one can put all frequencies to be zero 
since they provide only quantum corrections to the parti- 
tion function. The remaining part of the partition func- 
tion is determined solely by a potential energy and does 
not depend on velocities. The scenario becomes dif- 
ferent for a position-dependent mass. In this case one 
can use the following general arguments. The parti- 
tion function Z , which is proportional to the phase vol- 
ume Ap Ax, acqu ires an additional positional dependence 
Ap ~ ^jTm{x) followin g fro m the momentum channel 
even in the thermal limit [33|. The free energy (— TlnZ) 
results in the interaction potential pS]) obtained by the 
rigorous derivation procedure. 

In a multidimensional case the kinetic energy is ex- 
pressed through the mass tensor rriy (i?) as 
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As shown in Ref. [33j , in the multidimensional case the 
interaction due to variable mass has the form 
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FIG. 1: Arrangements of particles with velocities u to calcu- 
late in Eq. (gSJ the mass (a) and (b) A/^ . 



The potential ([201) has a fluctuation origin and is medi- 
ated by fast fluctuations of velocity. In terms of coordi- 
nates, when the form p9|) is diagonalized, det m becomes 
a product of principal values and the interaction (|20p is 
reduced to a sum of terms related to principal coordi- 
nates. 



VI. INTERACTION MEDIATED BY THERMAL 
FLUCTUATIONS OF PARTICLE VELOCITIES 

To calculate the effective fluctuation potential even for 
system with a complicated dynamics, one has to find in 
the high frequency limit the mass tensor and to insert 
it into Eq. (|20p . When particles in a fluid perform an 
oscillatory motion with a high frequency, the fluid ve- 
locity obeys the Euler equation everywhere in the fluid 
excepting a thing layer close to the particle surfaces [13] . 
Hence for finding the mass tensor in Eqs. (fT9)l - ((20)) one 
has to solve the Euler equation with the zero boundary 
condition for a normal component of the fluid velocity. 
For this reason, the mass, corresponding to the high fre- 
quency limit of particle dynamics, can be called Euler 
mass. The effective particle masses depend on a fluid 
mass involved into the motion. The fluid mass depends 
on an interparticle distance and therefore the effective 
particle masses also depend on that distance. 

In the case of one particle of a radius a in a bulk fluid 
the Euler mass tensor has the form 13411 



(21) 



(20) 



where po is a mass density of the particle and p is a fluid 
density. The first term in Eq. (PT|) is related to the self 
mass of the particles and the second one is associated 
with a fluid motion. 



A. Two particles in an infinite fluid 

Now we calculate the Euler mass tensor for two iden- 
tical particle in a fluid. The fluid velocity, normal to a 
particle surface, should equal a normal particle velocity 
on its surface. If the particle velocities are ui and U2 there 



FIG. 2: Arrangements of particles with velocities u to calcu- 
late in Eq. (|23|l the mass (a) irix and (b) m^. 



are four independent quadratic combinations {u\ + u^), 

U1U2, {Rui){Ru2), and [[Rui)'^ + {Ru2)'^]- Let us make 
the velocity transformation 



V = 
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In terms of new velocities the kinetic energy can be writ- 
ten in the form 



m,{R)vl 
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where i = 1, 2, 3 correspond to x, y, z respectively. In the 
expression (|23p there are four independent masses only 
since = My and nix — "niy- One can separate self 
masses of the particles 47ra'^/Oo/3 from fluid ones using 
the form 
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According to Eq. (pij) . G(oo) = g{oo) = 1. The easiest 
way to calculate the masses is to use the method illus- 
trated in Figs. [1] and [21 where the particle velocities are 
shown by the arrows. The kinetic energy of the fluid 
in Figs. [TJa) and[Hb) are {2TTa^ /3)Gx.ziR/a)u^ respec- 
tively. Analogously, the fluid kinetic energies in Figs. a) 
andlD^b) are /3)gx,z{R/ aju^ . The interaction po- 

tential (I20I) takes the form 
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The function I{R) tends to zero at i? — > 00. 



B. One particle near an infinite wall 



FIG. 3: The particle near two perpendicular walls. The three 
image particles, shown by the dashed curves, are added to 
consider the whole space with zero normal velocities of the 
fluid at the walls. This particles arrangement contribute to 
the function /i in Eq. (pf]!. 



center-to-plane (the plane is at z = 0) distance is h. The 
boundary condition at the particle surface is an equality 
of normal components of the fluid and the particle veloc- 
ities. The normal velocity component of the fluid at the 
flat surface is zero. Obviously, the total kinetic energy 
is one half of that calculated in the previous subsection 
corresponding to Fig.[IJa) (the x and y components) and 
Fig. [lib) which are related to the zero normal velocity of 
the fluid at = 0. 

It is easy now to write down the interaction potential 
using the results of the previous subsection. It reads 
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(26) 

where the functions Gx and g^ are the same as calculated 
in the previous subsection. 



One particle near two infinite perpendicular 
walls 



This situation is shown in Figs. [3]- [H The fluid is re- 
stricted by the conditions y < and z > where the 
particle is placed. The other image particles, figured out 
by dashed curves, serve to obey the boundary conditions 
in the planes z — and y = of zero normal velocities 
of the fluid. If to consider the whole space with incorpo- 
rated image particles the total kinetic energies are 
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with i = 1, 2, 3 for Figs. [3] - [5] respectively where the ve- 
locities of all particles are u. Analogously to the previous 
cases, the interaction potential is 



T^^^^2po + pMV^^. 
2 2po + p 



(28) 



We consider here one particle of the radius a placed 
in a half-space of fluid filled out the volume z > 0. The 



The boundary conditions are of the same type as in 
Sec. lvra 
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FIG. 4: The particles arrangement contribute to the function 
/2 in Eq. 



FIG. 5: The particles arrangement contribute to the function 
/3 in Eq. (HZJ. 



VII. NUMERICAL METHOD TO CALCULATE 
THE INTERACTION POTENTIAL 

In our case the fluid velocity is v = V(p, where the 
potential (p obeys the Laplace equation V^i^ = with 
the boundary conditions specified in Sec. IVII 



n{u — V(^) „ = 0. 



(29) 



In Eq. ([29)1 S is the total area restricted the fluid includ- 
ing particles, with velocities m, and walls with zero veloc- 
ities. The unit vector n is perpendicular to the surfaces. 
The kinetic energy of the fluid is / (firpiP' /2. Therefore, 
the total kinetic energy of the system can be written as 
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where summation is performed over all particles with the 
areas S^. 

The numerical method proposed to calculate Euler 
masses is based on an iteration procedure. Below we 
discuss this method for two particles in an infinite fluid. 
The zero approximation is 

if^{r)^,l:{r~R/2)+i!{r + R/2), (31) 

where the harmonic function ij: satisfies the equation 
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at the exterior of the sphere of the radius a and with the 
center at r = 0. At the sphere surface nWip{r) = ucosO 



and 6 is the angle between n and u. The integration in 
Eq. (|32p is extended over the sphere surface and Green's 
function has the form 
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and cos 7 



where |r — 7^'| = (r' + a' 
cos 9 cos 9' + sin 9 sin 9' cos( 

The potential fo(j) in the zero approximation (j3ip 
provides the correct boundary conditions at the parti- 
cle surfaces S in the limit R ^ 00. At a finite R the 
boundary condition (j29p is not satisfied by the zero ap- 
proximation ()31|) which requires a modification. One can 
construct the iteration scheme if ^ (po + cpi + ip2 + ■■■ by 
means of the recursion relation 

<p„+i (r) = - / -^G{f- r' - RI2) (34) 

J\r'-R/2\=a 47r 

n'{u5no - V(/3„(r")) 
R/2,r' + R/2) 



\r'+R/2\=a 
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n'{u5no ~ V(pn{r')), 



where n = 0, 1, 2, . 



At a large R there is a fast conver- 
gence since (fn+i ~ 'PnO-fR. At each iteration step the 
boundary conditions becomes more and more exact with 
respect the parameter a/R. 

To calculate = Gy and Gz in Eq. for Euler 
masses one has to apply the scheme (p4|) to the situations 
shown in Fig. [H^a) and in Fig. [Hb). Analogously, = 
gy and gz in Eq. ((26)l are associated with Fig. [U^a) and 
Fig.lHb). 

By means of the functions G^ and gz, numerically 
found by the above method, one can also construct the 
interaction (|26p of a particle and a wall. 

The same method is applicable to calculations of the 
potential ([28]) in Sec. IVII C for one particle near two per- 
pendicular walls. Instead of two particles one should take 
four ones when three of them play a role of images and 
are shown by the dashed curves in Figs. [3] - [S] We do 
not describe an obvious modification of Eq. ([5^ for that 
case. 



VIII. RESULTS 

In this chapter we discuss the three different contribu- 
tions to interaction of colloidal particles, listed in SecUl 
in various geometries: (i) two particles in a bulk fluid, 
(ii) one particle over a flat surface, and (iii) one particle 
near two perpendicular planes. 

The interaction Ucom, which is mediated by compres- 
sion fluctuations of the fluid, plays an outstanding role 
since it is relatively large, Ucom/T ~ (0.3 — 1.0), as fol- 
lows from Figs. [S] and [51 On the other hand, the interac- 
tion Ucom is strongly reduced by a finite Debye screening 
length Xo- Therefore, this interaction can be observed in 
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FIG. 6: The numerically calculated attraction potentials for 
two particles in an infinite fluid, p/po ~ 1. 
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FIG. 7: The numerically calculated attraction potentials for 
one particle and the infinite flat wall, p/po ~ 1- 



1.1 1.2 1.3 h/a 1.4 

I ■ ■ — 




-1.6 



FIG. 8: The thick curve is the same as in Fig. [T] The thin 
curve represents the interaction (|15p mediated by thermal 
compression fluctuations of the fluid. 




electrolytes with a large Xd- In addition to that, surfaces 
of interacting objects should not be strongly charged to 
prevent domination of the Coulomb repulsion over Ucom- 

In a conventional electrolytes, normally used in exper- 
iments, only Uydw (frequencies larger than plasma one) 
and / (non-compressive fluctuations) survive on distances 
larger than Xd where the Coulomb repulsion of charged 
particles is screened. 

Fig. [6] relates to the case of two particles in a bulk 
fluid. Sec. IVII A. The interaction potential / (in the units 
of temperature T) is plotted by the thin solid curve. The 
conventional van der Waals interaction ^ is indicated 
by the dashed curve. The resulting potential is shown by 
the thick curve. It can be seen that Ucom substantially 
exceeds the above interaction potentials. 

Figs. [7] and [5] correspond to the case of one particle 
near the wall. Sec. IVII B. As clear from Fig. [U Ucom is 
not small and is of the order of T. 

In Fig. [5] the potential / relates to one particle near 
two perpendicular walls. 

When the intersurface distance is small the potential 
/ tends to a constant whereas the van der Waals inter- 
action diverges as known. The results shown in Figs. [H] 
- [S] correspond to p/po = 1. Under increase of this pa- 
rameter the attraction becomes more pronounced. For 



FIG. 9: The attractive interaction / in the units of T for one 
particle and two perpendicular walls, p/po ~ 1- 

example, for p/ po = 13.6, related to colloidal particles in 
mercury, the attraction potential becomes roughly more 
than two times larger in amplitude. 

IX. DISCUSSIONS 

Is it possible to experimentally observe an attraction of 
like-charged colloidal particles separated by micron dis- 
tance? Let us analyze different physical situations. 

A. No Coulomb effects 

A lack of Coulomb effects means use of an electrolyte 
with a large Debye length Xd and zero charges on all 
surfaces. In this case DUVO repulsion is absent and the 
interaction mediated by compression fluctuation, Ucom, 
substantially dominates the conventional van der Waals 
part Uvdw and one mediated by particle velocities, /. 
As follows from Fig. [SJ the interaction energy Ucom of a 
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particle of the radius 1 fx and a flat plane, when surface- 
to-surface distance is 0.15 fx, is IT. This attraction can 
be experimentally observed since it is relatively large. 

B. Conventional particles in an electrolyte 

For typical electrolytes and polystyrene particles [IHl 
DLVO repulsion is of a short length and the interaction 
Ucom is well suppressed. In this situation only u^dw and 
/ survive. But the van der Waals interaction Uvdw is 
weak and only / has chances to result in an observable 
attraction. 

The peculiarity of the interaction / is that a particle is 
attracted stronger to a surface which better geometrically 
adjusts its shape. The flat surface in Fig. [7] adjusts the 
particle better than the counter-curved neighbored par- 
ticle in Fig. [51 Also the geometry of two walls in Fig. [5] 
adjust the particle shape better than one wall in Fig. [T] 
An attraction of a particle to a surface of the same type 
of curvature is most stronger. 

The feature of the interaction / is that it is of the order 
of 0.05T (a particle near a flat wall) and of the order of 
O.IT (a particle near two perpendicular walls) when the 
intersurface distance is of the order of 0.2a ^ 2000A. At 
a larger intersurface distances the attraction / exceeds 
the van der Waals part. 

It is promising to take a surface which adjusts the 
spherical particle better than two perpendicular planes. 
For example, it can be a cylinder with the particle in- 



side. In this case an interaction is a few or even ten 
times larger than O.IT (two perpendicular planes). The 
minimum value of an interaction potential O.IT, when T 
is room temperature, is approximately a border of ex- 
perimental resolution. So in a conventional electrolyte 
the attraction a micron sized particle to two perpendic- 
ular planes or to an interior surface of a cylinder can be 
experimentally observed. 



X. CONCLUSIONS 

For like charged colloidal particles two mechanisms of 
attraction between them survive when the interparticle 
distance is larger than the Debye screening length. One 
of them is the conventional van der Waals attraction and 
the second one is the attraction mechanism mediated by 
thermal fluctuations of particle position. The latter is 
related to an effective variable mass (Euler mass) of the 
particles produced by a fluid motion involved. The most 
stronger attraction potential (up to the value of tempera- 
ture T) corresponds to a situation of uncharged particles 
and a relatively large Debye screening length. In this 
case the third mechanism of attraction enters the game. 
It is mediated by thermal fluctuations of fluid density. 
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